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Finite-size effects in the mean-field Ising spin glass and the mean-field three-state Potts glass are 
investigated by Monte Carlo simulations. In the thermodynamic limit, each model is known to 
exhibit a continuous phase transition into the ordered state with a full and a one-step replica- 
symmetry breaking (RSB), respectively. In the Ising case, Binder parameter g calculated for various 
finite sizes remains positive at any temperature and crosses at the transition point, while in the Potts 
case g develops a negative dip without showing a crossing in the g > region. By contrast, non-self 
averaging parameters always remain positive and show a clear crossing at the transition temperature 
in both cases. Our finding suggests that care should be taken in interpreting the numerical data of 
the Binder parameter, particularly when the system exhibits a one-step-like RSB. 

PACS numbers; 05.50.+q, 75.50.Lk 



I. INTRODUCTION 

The concept of replica-symmetry breaking (RSB) ||] 
gives us new insight into the character of the ordered 
state of complex systems such as spin glasses (SG) |^ and 
real structural glasses Systems exhibiting the RSB 
can roughly be divided into two categories depending on 
their breaking patterns: One is a full or hierarchical RSB 
and the other is a one-step RSB. In both cases, there 
are many different equilibrium states unrelated by global 
symmetry of the Hamiltonian, and an overlap between 
these states plays an important role in describing the 
ordered state. 

In the case of one-step RSB [Q , the overlap q takes only 
two values in the thermodynamic limit, namely, either a 
self-overlap equal to the Edwards- Anderson order param- 
eter, q — qEA, or a non-self-overlap usually equal to zero, 
q = 0. The overlap distribution function P{q) consists 
of two distinct delta-function peaks, one at g = gsA and 
the other at g = 0. One-step RSB transitions could be 
either continuous or first-order, either with or without a 
finite discontinuity in ^ea at the transition. Examples of 
the first-order one-step RSB transition may be the mean- 
field p-spin glass with p > 2, the random energy model, 
and the mean-field p-state Potts glass with p > 4, while 
those of the continuous one-step RSB transition may be 
the mean- field p-state Potts glass with 2.8 < p < 4. 

In the case of the full RSB, by contrast, possible values 
of the overlap are distributed continuously in a certain 
range, and the states are organized in a hierarchical man- 
ner. The overlap distribution function has a continuous 
plateau at q < q^A in addition to the delta-function peak 
at q = qEA- Well-known example of this category is the 
standard mean-field Ising SG, namely, the Sherrington- 
Kirkpatrick (SK) model. In some special cases, the ad- 
mixture of the above twos, where the overlap distribu- 
tion function has a continuous plateau together with the 



delta-function peak at q = (and the one at q — qEA), is 
also possible. An example of this may be the mean-field 
p-state Potts glass with 2 < p < 2.8. 

Recent interest in SG studies has been focused largely 
on the validity of applying the RSB idea established in 
some mean-field models to more realistic short-range SG 
models. In almost all such studies, the three-dimensional 
(3D) Ising SG model has been employed. Indeed, some 
researchers have claimed that the full or hierarchical RSB 
as observed in the SK model Q is also realized in realistic 
3D SG while other researchers have claimed, based on 
an alternative droplet picture that the ordered state 
of realistic 3D SG is unique up to global symmetry of 
the Hamiltonian, without showing RSB of any kind 
Thus, intensive debate has continued between these two 
scenarios as to the true nature of the SG ordered state 
of 3D short-range systems. 

Meanwhile, the one-step RSB has been discussed 
mainly with interest in its close connection to struc- 
tural glasses rather than SG magnets [§ 10|. Recently, 
however, one-step RSB features have been found unex- 
pectedly by the present authors in the chiral-glass state 
of a 3D Heisenberg SG [0. According to the chiral- 
ity mechanism of experimental SG transitions based on 
the spin-chirality decoupling-recoupling scenario fl^ , the 
SG ordered state and the SG phase transition of real 
Heisenberg-like SG magnets possessing weak but nonzero 
magnetic anisotropy are governed by the chirality order- 
ing of the fully isotropic system which is "revealed" by 
the weak magnetic anisotropy, not by the spin ordering 
which has been "separated" in the fully isotropic case 
from the chirality ordering. Then, the observation of 
Ref. g means that the SG ordered state of most of real 
SG magnets should also exhibit such one-step RSB-like 
features. Note that such a picture of the SG ordered 
state contrasts with the standard pictures discussed so 
far, either the droplet picture without the RSB or the 
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SK picture with the full RSB. 

Under such circumstances, further studies of the na- 
ture of the possible RSB in 3D short-range SG models 
are clearly required. Since we are usually forced to em- 
ploy numerical simulations to investigate 3D short-range 
models, and since numerical simulations are often ham- 
pered by severe finite-size effects, we feel it worthwhile to 
further clarify by numerical simulations the finite-size ef- 
fects in some mean-field models which are exactly known 
to exhibit RSB transitions in the thermodynamic limit. 
In particular, the question how the one-step and full 
RSB transitions look like in finite-size simulations is of 
both fundamental and practical interest. Such informa- 
tion would be of much help as a reference in interpreting 
the numerical data obtained for finite-dimensional short- 
range SG models. 

In the present paper, we choose two mean-field SG 
models exactly known to exhibit a continuous (second- 
order) phase transition in the thermodynamic limit: One 
is the SK model which shows the full RSB, and the other 
is the mean-field three-state {p = 3) Potts-glass model 
which shows the one-step RSB. We calculate by Monte 
Carlo simulations several quantities which have widely 
been used in identifying the phase transition, including 
the spin-glass order parameter and the Binder parame- 
ter, together with the quantities recently introduced to 
represent the non-self-averaging character of the ordered 
state. By carefully examining the size dependence of 
these quantities, comparison is made between the two 
types of RSB. Our results have revealed that the Binder 
parameter of the one-step RSB system shows the be- 
havior very different from the standard behavior, giving 
warning about the interpretation of the numerical data 
for relevant short-range systems. 



II. MODELS 

The mean-field p-state Potts-glass model is defined by 
the Hamiltonian, 



n = 



N 

i<j 



(1) 



where Ui denotes a Potts-spin variable at the i-th site 
which takes p distinct states, and N is the total number 
of Potts spins. The exchange interaction is an inde- 
pendent random Gaussian variable with zero mean and 
variance J'^/N. The model with p = 2 is equivalent to 
the SK model. In the present study, we focus our atten- 
tion on the standard SK model corresponding to p — 2 
and the three-state Potts-glass model corresponding to 
p = 3. Although the thermodynamic properties of an 
infinite system have been rather well understood by the 
calculation based on a replica technique its 
finite-size properties have been much less understood. 



It is convenient to use an equivalent simplex spin rep- 
resentation where the Potts spin Ui is written in terms 
of a p — 1 dimensional unit vector Si, which satisfies 
Si • Si 



p-i 



N 



H = -{p-l)Y^J,,S,-S,. 



(2) 



In the particular case of p — 2, Si simply reduces to the 
one-component Ising variable Si = ±1, and the Hamilto- 
nian (2) is equivalent to the standard SK Hamiltonian. 

In terms of the simplex spin >5'f(l</i<p— 1), the 
parameter q may be defined by 



(3) 



where g'"' denotes an overlap tensor between two replicas 
1 and 2, 
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TV ^ 
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The Binder parameter is then given by 

2 m] 



{p-lf 



(4) 



(5) 



where (■ • •) denotes the thermal average and [■ • ■] denotes 
the average over the quenched randomness {Jy}. The 
Binder parameter is normalized so as to vanish above the 
transition temperature Tg in the thermodynamic limit. 
Recall that at T > Tg each component q^^i, should be- 
have as an independent Gaussian variable. Below Tg, g 
is normalized to give unity in the thermodynamic limit 
for the nondegenerate ordered state where P{q) has only 
trivial peak at g = (7ea- Of course, this is not the case 
for the SG models showing the RSB including the present 
mean-field SG models, for which g takes nontrivial val- 
ues different from unity even in the thermodynamic limit. 
Hence, at least in the case where a continuous phase tran- 
sition occurs into the trivial ordered state, g for various 
finite sizes is expected to cross at T = Tg. Indeed, this 
aspect has widely been used for locating the transition 
temperature from the numerical data for finite systems. 



III. MONTE CARLO RESULTS 

We perform MC simulations based on a version of the 
extended ensemble method, called the exchange method 
Jisf . As in other SG models, an extremely slow relax- 
ation becomes a serious problem of MC simulations in 
the present mean-field SG models. Such difficulty could 
partly be overcome by using the exchange method, which 
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has turned out to be quite efficient in thermahzing vari- 
ous hardly relaxing systems. The method enables us to 
study larger sizes and/or lower temperatures than those 
attained previously. Our MC simulations have been per- 
formed up to = 512 at T/J = 0.25 for the SK model, 
and TV = 256 at T/J = 0.4 for the mean-field p = 3 Potts- 
glass model, where Tg/ J = 1 in both models. Sample av- 
erages are taken over 200 — 1792 independent bond real- 
izations depending on the size N . We note that the min- 
imum temperature studied here are considerably lower 
than the previous ones; e.g., N = 512 at T/J = 0.75 ||l^] 
for the SK model and N = 120 at T/J = 0.98 for the 
mean- field p = 3 Potts-glass model [|7| . 

The temperature and size dependence of the calculated 
Binder parameter g is shown in Figs, [l] and § for the 
SK and the p — 3 Potts-glass models, respectively. As 
is evident from these figures, the Binder parameters of 
the two mean-field models show considerably different 
behaviors from each other. 

In the SK model, as shown in Fig. |^, a clear crossing of 
g is observed at T = Tg, which looks similar to the ones 
seen in the standard continuous transitions. In fact, the 
behavior of g found here also resembles the ones observed 
in the short-range Ising SG models in 3D jlj,|l^ and in 
4D |2^ , though the crossing tendency is less pronounced 
in 3D than in 4D. As mentioned, g of the SK model takes 
a nontrivial value below Tg even in the thermodynamic 
limit due to its RSB. We show in Fig. |l] the behavior 
of g{T, oo) evaluated in the replica formalism by numer- 
ically solving the Parisi equation [1(a)]. Note that, as 
the temperature approaches Tg from below, the limiting 
value g{T~ , oo) goes to unity as in the case of ordinary 
continuous phase transitions. Hence, with increasing N , 
g{T~ , N) just below Tg is expected to approach unity 
from below while g{T^ , N) just above Tg approaches zero 
from above, which entails a crossing of g at T = Tg. With 
lowering the temperature, g{T, oo) first decreases, reach- 
ing a minimum around T/J — 0.5, and increases again 
tending to unity at T = 0. Here note that, for any model 
with nondegenerate ground state, P{q) becomes trivial 
at T = irrespective of the occurrence of RSB, and g 
tends to unity. As can be seen in Fig. 0, the present MC 
results for finite N gradually approach the g{T, oo) curve 
of an infinite system. 

In the mean-field p — 3 Potts glass, as shown in Fig. ||, 
no crossing of g is observed at T = Tg , at least of the 
type as observed in the SK model. Instead, unlike the 
case of the SK model, a shallow negative dip develops 
above Tg for larger N which becomes deeper as the sys- 
tem gets larger. Although the existence of a negative dip 
was not reported in the previous numerical works [ pT| , we 
note that a negative dip appears only for larger N which 
accounts for the absence of a negative dip in the previ- 
ous data. Perhaps, on looking at Fig. |^, one would hardly 
imagine that there occurs a continuous phase transition 
at T/J = 1: Nevertheless, the occurrence of a continuous 



transition at T/J = 1 is an exactly established property 
of the model. We also note that, while the appearance 
of a growing negative dip in the Binder parameter is of- 
ten related to the occurrence of a first-order transition 
p2[ , this is not always the case: Here, the transition is 
established to be continuous. 

It might be instructive to examine here the behavior 
of g in the thermodynamic limit. As the temperature 
approaches Tg from below, giT~ , oo) tends to a negative 
value, —1 in the present case. Such a negative value of 
g{T~ , oo) is in sharp contrast to the system showing the 
full RSB where g{T~ , oo) = 1. Indeed, this negativity is 
closely related to the occurrence of the one-step RSB in 
the model [||. 

Then, one expects that the negative dip of g{T, N) 
observed in Fig. further deepens with increasing N, 
and eventually approaches —1 from, above at T — T^ , in 
sharp contrast to the SK case where g(T<r , N) approaches 
1 from below. Therefore, the crossing of g in the g > 
region as observed in the SK model hardly occurs in the 
p = 3 Potts-glass model. Rather, if one considers the 
fact that g(T, TV) above Tg is negative for moderately 
large N approaching zero from below, the crossing of g is 
expected to occur in the g < region, not in the 5 > 
region as in the case of the SK model. The data of Fig. | 
are certainly consistent with such a behavior. Anyway, 
our present result of the mean-field p = 3 Potts glass 
has revealed that the data of the Binder parameter has 
to be interpreted with special care particularly when the 
ordered state has one-step RSB features. 

Next, we study the so-called Guerra parameter which 
was originally introduced to detect the RSB transition 



G{T,N) 



(6) 



Since the numerator represents a sample-to-sample fluc- 
tuation of the overlap, non-vanishing of G means a lack 
of self-averaging so long as the denominator remains 
nonzero. In the mean-field SG models studied here, their 
RSB indeed gives rise to the lack of self-averaging, i.e., 
the occurrence of a non-trivial probability distribution 
of the overlap over quenched disorder. It has been rig- 
orously proven, without using the replica trick, that in 
the SG phase of the SK model the G parameter in the 
thermodynamic limit is equal to 1 /3 independent of the 
temperature [ p5| . Meanwhile, it has been pointed out in 
Rcf. that, even when P{q) is trivial and the ordered 
state is self-averaging, the G parameter can still take a 
nonzero value due to the possible vanishing of the de- 
nominator, leading to a crossing at Tg |26|. Hence, the 
crossing of G does not necessarily mean the lack of the 
self-averaging, although it can still be used as an indi- 
cator of a phase transition. As an indicator of the non- 
self-averageness in the ordered state, one may use the A 
parameter defined by PTl 
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cal dimension vmf, £, ^ \T — Tg\ '^'^^ . Then, the squared 
order parameter can be written as 



We calculate these two parameters, G and A, both for 
the SK and the mean-field p = 3 Potts-glass models. The 
temperature and size dependence of the G and A param- 
eters of the SK model is shown in Figs. || and |[ respec- 
tively. Although the error bars are still large, both G and 
A show a clear crossing at Tg, remaining positive at any 
temperature. As expected, with increasing TV, the G pa- 
rameter approaches 1/3 independent of T below Tg. By 
contrast, the A parameter for various sizes merge into 
a curve below Tg, which clearly stays nonzero indicat- 
ing the non-self-averageness of the ordered state. Here 
it should be noticed that, just at the transition point, 
the non-self-averageness is expected to occur in any ran- 
dom system, even includin g th e ones without showing the 
RSB in the ordered state |28y29| ] . Hence, in the type of 
random systems which do not show the RSB in the or- 
dered state, A(T, oo) stays nonzero only just at T = Tg 
and vanishes on both sides of Tg. By contrast, in the 
present SK model, A(T, oo) should stay finite even below 
Tg due to its RSB, which explains the observed merg- 
ing behavior seen in Fig. ^ at T < Tg. As can be seen 
from Fig. ^, on further lowering the temperature toward 
T — 0, A{T, N) tends to vanish in contrast to the behav- 
ior of G{T,N). This aspect is consistent with the fact 
that at T = the overlap distribution becomes trivial 
and the self-averageness is recovered irrespective of the 
occurrence of RSB. 

The G and A of the mean-field p — 3 Potts glass are 
presented in Figs. ^ and ||, respectively. Unhke the case 
of the Binder parameter g shown in Fig. ||, the G and 
A parameters remain positive at any T and show a clear 
crossing at T = Tg: They behave more like the Binder 
parameter of standard systems, e.g., like the one shown in 
Fig. 1^. In fact, the behaviors of the G and A parameters 
shown in Figs, [s] and pi are similar to those of the SK 
model shown in Figs, g and ^, suggesting that G and 
A are less sensitive to the kind of breaking pattern of 
replica symmetry. Hence, one could use the G and A 
parameters to identify the SG transition based on the 
standard crossing method even for systems showing a 
one-step RSB. 

Once the transition temperature is established, the 
next task would be to determine critical exponents. Here 
we wish to examine a finite-size scaling hypothesis con- 
cerning the SG order parameter for the present mean- 
field models. Similar analysis has widely been used 
for extracting the critical exponents from the numeri- 
cal data. According to Ref. , finite-size scaling of the 
mean-field models can be derived by assuming that the 
"coherence number" behaves as where du is the up- 
per critical dimension of the corresponding short-range 
model, while the "coherence length" ^ diverges at T = Tg 
with the correlation-length exponent at the upper criti- 



[{q^)] - |T - Tg|2/3MF j(jV|T - Tgl''"''"''), 

~ N-'^P^F/d^i^mF f'(^]S[\T — Tgl''"""''), (8) 

where /3mf = 1 is the mean-field order-parameter expo- 
nent whereas / and /' are the scaling functions. Not- 
ing the fact that the upper critical dimension of the SG 
models is = 6 and the correlation-length exponent at 
d = d„ = 6 is equal to j'mf=1/2, it follows 

^ N-^/'^fWT - Tg\N^'^). (9) 

The resulting finite-size scaling plots are shown in Figs. 
and H for the SK and the p = 3 Potts-glass models, re- 
spectively. In both models, the scaling of the form (||) 
turns out to work fairly well both below and above Tg 
as far as the temperature is sufficiently close to Tg. We 
note that a similar finite-size-scaling analysis has already 
been reported for the SK model just at T„ |3l| and for 
the p = 3 Potts-glass model above Tg [Trf . In particu- 
lar, the scaling turns out to be reasonably good even for 
the p = 3 Potts glass where the Binder parameter does 
not exhibit a clear crossing in the range of sizes studied. 
This implies that the standard finite-size scaling analysis 
of the order parameter could still be useful even in RSB 
systems including the one-step RSB systems. 

IV. DISCUSSION AND REMARKS 

In this section, with our present results for the mean- 
field models in mind, we wish to comment on the possible 
RSB in some short-range SG models. 

As mentioned, one-step RSB-like features were re- 
cently observed in the chiral-glass state of the 3D short- 
range Heisenberg SG There, the Binder parameter 
for the chirality, the order parameter of the chiral-glass 
transition, did not cross in the g > region and de- 
veloped a negative dip which deepened with the system 
size. Instead, a crossing of g was observed in the g < 
region close to the negative dip (see Fig. 1 of Ref. 
). Meanwhile, the G parameter always remained posi- 
tive and showed a clear crossing at T = Tg (see Fig. 3 
of Ref. ). All these features are similar to the ones 
observed here in the mean-field p = 3 Potts glass, sug- 
gesting that the chiral-glass state of the 3D Heisenberg 
SG has a one-step RSB-like character |Q . 

Other obvious interest is the nature of the possible 
phase transition of the short-range three-state {p = 3) 
Potts-glass model in 3D. It is widely believed that there is 
no finite-temperature phase transition in 3D p = 3 Potts 
glass which were investigated by MC simulations ]3^-|3^] 
and other numerical methods ||3^. In particular, MC re- 
sults of Refs. ||3^,^ revealed that the Binder parameter 
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decreased monotonically with system size without show- 
ing a crossing, which was taken as an evidence of the 
absence of a finite-temperature transition. However, the 
behavior of g observed in Refs. [^,^ was not dissimilar 
to the one observed here in the mean-field p = 3 Potts 
glass, and we feel that the possibility of the occurrence 
of a one-step RSB-like transition at finite Tg still cannot 
be ruled out. 

Recently, short-range p-spin glass models whose mean- 
field versions have been known to show the one-step RSB 
were studied by MC simulations p7|,p8[ . For example, 
according to the calculation of Ref. |j8[ for the 4D 3-spin 
model, the Binder parameter did not exhibit a crossing of 
the standard type, while the G and A parameters showed 
a clear crossing at T = Tg > 0, strongly suggesting the 
occurrence of a finite-temperature transition. Thus, from 
out present study, the possible occurrence of a one-step 
RSB transition at T = Tg > is suspected. Meanwhile, 
a closer inspection reveals that a negative dip observed in 
g becomes shallower with increasing system size in 
contrast to the case of the mean-field p = 3 Potts glass 
studied here. Further studies seems to be required to 
clarify the nature of the RSB in the short-range p-spin 
glass. 

In conclusion, we have investigated by MC simulations 
the finite-size effects of the two mean-field SG models 
whose replica-symmetry-breaking properties in the ther- 
modynamic limit are well established. In the mean-field 
Ising spin glass (the SK model), the Binder parameter 
g of various sizes always remains positive and crosses at 
T = Tg, while in the mean-field three-state Potts glass, 
it develops a negative dip which deepens as the system 
size increases, without a crossing in the g > region as 
observed in the SK model. Instead, a crossing of g occurs 
in the g < region near the negative dip. Such difference 
in the behaviors of g reflects the different types of asso- 
ciated RSB of the two models, i.e., full versus one-step 
RSB. By contrast, the Guerra parameter G as well as 
the non-self-averaging parameter A always remain posi- 
tive and show a crossing of the standard type at T = Tg 
for both the SK and p — 3 Potts-glass models. We have 
also discussed implications of the present results to the 
possible interpretation of the numerical results for some 
short-range SG models. 
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FIG. 1. Temperature and size dependence of the Binder 
parameter of the SK model. The bulk transition temperature 
is located at T/J = 1. The broken line represents the RSB 
solution derived by solving the Parisi equation. 
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FIG. 2. Temperature and size dependence of the Binder 
parameter of the mean-field p — 3 Potts glass. The bulk 
transition temperature is located at T/J = 1. 
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FIG. 3. Temperature and size dependence of the G param- 
eter, defined by Eq. (^) , of the SK model. The bulk transition 
temperature is located at T/J — 1. The broken straight line 
represents the line G — 1/3. 
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FIG. 4. Temperature and size dependence of the A param- 
eter, defined by Eq. (^ , of the SK model. The bulk transition 
temperature is located at T/J = 1. 
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FIG. 5. Temperature and size dependence of G parameter, FIG. 8. Finite-size scaling plot of the squared order pa- 

defined by Eq. (^), of the mean-field p = 3 Potts glass. The rameter of the p = 3 mean-field Potts glass with the scaling 
bulk transition temperature is located at T/J = 1. form Eq. with Tg/J = 1. 




FIG. 6. Temperature and size dependence of A parameter, 
defined by Eq. (^), of the mean-field p = 3 Potts glass. The 
bulk transition temperature is located at T/J — 1. 
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FIG. 7. Finite-size scaling plot of the squared order pa- 
rameter of the SK model with the scaling form Eq. with 
TJJ=1. 
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